Electron-Phonon Correlation Effects in Molecular Transistors 
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The interplay of electron-electron and electron-phonon interactions is studied analytically in the 
Kondo regime. A Holstein electron-phonon coupling is shown to produce a weakening of the gate 
voltage dependence of the Kondo temperature and may explain the observed anomalies in some of 
these devices. A molecular center-of-mass mode opens a new channel for charge and spin fluctuations 
and in the antiadabatic limit the latter are described by an asymmetric two-channel Kondo model. 
Below the Kondo temperature the system develops a dynamical Jahn- Teller distortion and a low 
energy peak emerges in the phonon spectral density that could be observed in Raman microscopy 
experiments. 

PACS numbers; 72.15.Qm, 73.22.-f 



I. INTRODUCTION 

One of the most exciting aspects in the physics of 
molecular transistors is the interplay between electron- 
electron and electron-phonon (e-ph) interactions/i 
Molecules attached to leads act like small quantum dots 
with broadened but well defined electronic levels and 
large charging energies. Electron-electron interactions 
lead to the well known Coulomb blockade effect^'^ and 
the strongly correlated Kondo physicsj^i^i^ When the 
coupling of electronic degrees of freedom and molecu- 
lar vibrational modes is strong, a rich variety of be- 
havior emergesi^i^^^^^i^^^^iii^i^^i^^il^^ The 
presence of steps in the current-voltage characteristics of 
molecular transistors can be explained by the coupling 
to a well defined phononic mode. Other features, like the 
anomalous behavior of the Kondo temperature as a func- 
tion of the gate voltage, remain yet to be understood!^ 
Center-of-mass oscillations of the active component of a 
transistor and electron-electron interactions are expected 
to produce a shuttling effect in nanomechanical devicesi^ 
The interplay of such mode with strong correlations at 
low temperatures is quite a complicated problem of gen- 
eral interest and experimental relevance^ 

Here we consider a molecule in the Kondo regime with 
a Holstein coupling and a center-of-mass motion mode 
(CMM) that modulates asymmetrically the coupling be- 
tween the molecule and the leads (see Fig. [T]). The 
rest of the paper is organized as follows: in next section 
we present the model and obtain the Kondo Hamiltoni- 
ans describing the low energy physics, in section HI we 
first analyze the behavior of Kondo couplings and Kondo 
temperatures, then we present variational wave functions 
that allow to explore a wider range of parameters, in par- 
ticular situations in which the phonon frequencies are of 
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FIG. 1: Scheme of the molecular device studied in this work. 
A molecule is attached to contacts: L and R, the center-of- 
mass motion modulates the molecule-leads tunneling barriers. 



the order or smaller than the Kondo temperature. Sec- 
tion IV contains a study of the spectral density of the 
CMM phonon. Finally, last section summarizes results 
and conclusions. 



II. MODEL 

The model Hamiltonian is H — Hnj + He + Hme 
where the first two terms describe the isolated molecule 
and electrodes, respectively, and the last term describes 
their coupling. We have 

Hm = £dnd + U Ud^ridi ~ X{nd~ 1) (a -I- a^) (1) 
He = ^"'^ ^i/caCafca. (2) 



Hme = ^ Vaixb) (dl 



(3) 



'Deceased 



Here rid = n-d^ + ndi, rida = dtda, dj- creates an electron 
at the molecular orbital with energy Sd, U is the intralevel 
Coulomb repulsion and c^^.^. creates an electron in the k 
mode of lead a = R, L. The operators and 5^ create a 
phonon excitation on the Holstein mode with energy ujq 
and on the CMM with energy oji , respectively. The single 
particle energies are measured from the Fermi energy and 
H=l. 

From hereon we consider a symmetric molecule with 
identical R and L leads {Sak = £k) and take Vji{xb) = 
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Vo + gxb and VL{xb) = Vo - gxb, where Xb = 
(l/2Mwi)^/^(& + b^) is the molecular center-of-mass co- 
ordinate. For a molecule with inversion symmetry, Xb 
cannot couple linearly with the electron density. 

It is convenient to rewrite the Hamiltonian in terms of 
symmetric and antisymmetric modes defined by cska = 

(cRkcr + CLka)/V2 and CAkcr = {cRk<7 - CLka)/V2: 



Hme = VsJ2i'^Uska + cl^^d„) 

+ VAib + b^)Y,idlcAk. + c\,^d^) (4) 



k.i7 



with Vs = V^Vo and Va = ^/TJWujlg. This shows that 
the symmetric mode is directly coupled to the molecu- 
lar orbital while the antisymmetric coupling is phonon- 
assisted. 

In the Kondo limit we eliminate the charge fluctuations 
by means of a Schrieffer- Wolff transformation extending 
the analysis of Refs. ^f. 12,21]. We obtain that the spin 
fluctuations are described by an asymmetric two-channel 
Kondo Hamiltonian 
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Here S is the spin operator of the molecular orbital, a 
are the Pauli matrices, and 
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-Ed + mujo Sd + U + mojo 



(6) 



2Vi7„ 



-Ed + + mwo Ed + U + UJI+ mujo 



where 7„ 



Ed + X^/UJQ, 



U = U — 2A^/a;o, £d <0, and U > 0. To second order in 
Hme there are also inelastic contributions proportional 
to VsVa which couple the ground state of the isolated 
molecule to a molecular state with a CMM excitation. 
The energy level structure of the molecule is analogous to 
that of a multilevel quantum dot with level splitting lui , 
and for large oji these contributions can be neglected;^^ 
We will present below a more detailed calculation where 
this assumption is not required. 



III. SPIN SCREENING IN THE KONDO LIMIT 

The low energy properties of Hk are well known .-^ 
At low temperatures only one channel participates in the 
screening of the molecular spin, and the Fermi liquid pic- 
ture applies as T ^ 0. In the renormalization group 
language, the zero temperature fixed point corresponds 
to a diverging ratio of the largest to the smallest Kondo 
coupling. Only for Js = Ja both channels participate in 
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FIG. 2: a) The Kondo temperature plotted against ei; for 
different values of the e-ph coupling A. X/u>o = 1 (thin line), 
2.5 (dotted line), 3.5 (dashed Une), and 4.5 (thick line). Here 
Ts = TTpVs = 20meV, U ^GTs, and loq/U = 0.5. b) Kondo 
coupling Js against Cd for the parameters of a). 



an overscreening of the molecular spin. In our context, 
this would be a quite peculiar situation. In the moderate 
or weak e-ph coupling regime Va is smaller than Vs and 
the Kondo coupling Js is therefore larger than J a- The 
screening of the molecular spin is then dominated by the 
symmetric channel. In this case the Kondo temperature 
is given by Tk ~ De^^^P'^^, where I? is a high energy 
cutoff and p is the bare electronic density of states of the 
leads at the Fermi level. These results imply that, at 
zero temperature, the asymmetric channel is effectively 
decoupled from the molecule and the conductance is gov- 
erned by the symmetric channel. As a consequence the 
conductance is simply 2e^ /h. 

In what follows we show that the obtained Js may 
explain some experimentally observed anomalies on the 
gate voltage dependence of Tk In the absence of e- 
ph coupling (A — 0) Js is given by the usual expression 
Js{^ = 0) = —2VsU /[Ed{Ed -f U)] and increases as the 
gate voltage decreases le^l or (e^ -I- U). The experimen- 
tally observed dependence of the Kondo temperature on 
the gate voltage is, however, much weaker than what is 
inferred from this expression. In Ref. 6 the Kondo tem- 
perature remains almost independent of the gate voltage 
and presents a fast increase very close to the charge de- 
generacy points. Other systems also present a change in 
the Kondo temperature dependence on gate voltage at 
the borders of the Coulomb blockade valleyi^ 

These anomalies can be understood by the presence of 
e-ph coupling. Using the expression for Js from Eq. ([6]) 
we plotted the gate voltage dependence of the Kondo 
temperature for different values of the e-ph coupling (see 
Fig. [2]). The gate voltage dependence of Tk weakens 
monotonically with increasing e-ph interaction. This can 
be seen more easily making some simplifying approxima- 
tions on Eq. ^ for Ed ^ —U/2. Note that the relevant 
quantities here are Ed and U since the valence changes 
in the molecule occur for = and Ed = —U. One 
of the effects of the e-ph coupling is to decreases both 
\£d\ and J7 shrinking the region of valence stability. For 
large A/wq 3> 1, 7(^1) = 7™ is approximately Gaussian 
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peaked at m ~ m* = (A/wo)^ of width ^/m* and satisfies 
X]m7m — 1- t^i^ case, the series in Eq. © can be 
approximated with 



Js 
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-Ed + m*ujo [/ + + m*ujo 
1 1 



U + Ed 



(7) 
(8) 



Wc obtain in the large X/loq limit that Jg is given by 
the A = expression! This result, however, ceases to be 
valid close to the charge degeneracy points Ed 0,—U 
where the m — term in Eq. ^ diverges with an expo- 
nentially small Franck-Condon prefactor 70 = e"''^/'^"'^ . 
The smallness of this prefactor requires Ed or ^d + U to 
be exponentially small in order to make the m = term 
of order one [see Fig. [2Ib)]. As a consequence, Js be- 
haves as Js(A = 0) which is regular at the charge degen- 
eracy points and only very close to them in behaves as 
Js ^ -2V§if /[£d{Ed + U)], which is divergent. The lat- 
ter expression gives the behavior of J5 that can naively 
expected using the renormalizcd single particle energy Ed 
and Hubbard interaction U in the usual expression for the 
Kondo coupling in an Anderson model. The behavior of 
the Franck-Condon matrix elements 7^ is the source of 
the anomalous dependence of the Kondo temperature on 
the gate voltage. Related anomalies are obtained in the 
valence change of simpler non-interacting model*^ and in 
the incoherent non-linear transportj^ 

It is important to notice that to lowest order in the 
hybridization only the Holstein coupling renormalizes Js 
and produces the anomalous behavior. In what follows 
we show that in higher order the CMM contributes to 
increase the Kondo coupling an effect that may help ex- 
plain the high values of Tk observed in some devices. To 
do this we now return to the original description given by 
the full Hamiltonian H and build the ground state wave- 
function using a generalized Varma-Yafet ansatz.^** While 
this method is exact only in the large N limit, where the 
spin symmetry is extended to the SU{N) group, it gives 
the correct zero-temperature physics even in the spin 1/2 
{N = 2) case. Following ReiM- we propose for the ground 
state wavefunction: 



1$) = 60 \(bF) \Po) + J2 ^-f^ l'^-*^) 1^-'=) 



(9) 



k,a,v=A,S 



where \4)f) describes the Fermi sea with total energy 
Efs. \(t)i^k) = TjEa^Cz^fcff l^o), and \P^k) are 

phonon states including both the Holstein coupling and 
the CMM, and we considered the U ^ 00 limit. The 
coefficients bo and buk are chosen for the wavefunction to 
satisfy the Schrodinger equation H \ ^) — Eqs \ ^)- 

A simple procedure allows to determine the form 
of the phonon states and their coefficients. First wc 
note that (0^^ \H\ $) = Ecsbo \Po) and \H\ $) = 
Eosbuk \Puk}- This generates a system of coupled 
equations for the phonon states that, after eliminating 



&i/fe \Pi^k), can be cast in the form 

iSE + Ed)\Po) = [H°j + 2V^T{E) 

+ 2Vi{b + b^)r(E){b + b^)]\Po).{lQ) 

Here 5E — Eqs — Eps — £d is the energy gained due 
to the hybridization, E = —SE ~ Ed + Hm, Hm 
projection of Hm on the subspace with € electrons in the 
molecular orbital, and 
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(11) 



where we considered for simplicity a square half-filled 
conduction band of total width 2D. We will study first 
the case with A = 0. In this case the state \Po) de- 
scribes only the CMM degrees of freedom and can be 
expanded in the base of eigenstates of the uncoupled 
phonon |Po) = J2n'^n\n)- The following recursion re- 
lation is obtained for the coefficients: 



MnnCtn = Mn,n+2an+2 + Mn,n-2an-2, (12) 



where 



M„ „ = 5E + Ed- uJin - 



2y|r„ 



- 2y2[(n + i)r„+i + nr„_i]^ 

Mn^n+2 = 2Vi^(n + l){n + 2)Tn+i, 
Mn,n-2 = 2VWnin-l)Tn-i, 



(13a) 
(13b) 
(13c) 



with r„ = T{—SE + uJin). The phonon state \Po) can 
only include even or odd n states and therefore has a 
well defined parity. We can find approximate solutions 
of these equations by noting that in the Kondo limit 6E 
is small. We have 
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2D 



■In 



SE 



D 



and r„ 



1 
2D 



In 



LUlTl 



D + LOin 



that give for the even-|Po) subspace 

{SE + Ed~ 2y|ro - 2VlTi ~ E5)ao 



for n > 0, 
(14) 

(15) 



0, 



where S5 is a small correction that can be written as 
a continued fraction after eliminating all coefficients 
with n > 0. Equation ([T5]) gives a simple self-consistent 
condition for SE and in the weak coupling regime the 
usual solution is obtained 



SE 



where p — {2D) ^ and 



J's 



-De "■'s , 



2F| 



\Ed\ + 2v^ri + j:s' 



(16) 



(17) 



The energy gain SE is associated to the formation of the 
Kondo singlet and is of the order of the Kondo temper- 
ature. This result shows that, in the even-|Po) subspace. 
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the molecular spin is Kondo screened by the symmetric 
channel, and the antisymmetric one only contributes to 
a renormalization of the molecular energy. This is the 
ground state for Vs > Va- In the opposite situation the 
ground state belongs to the odd-|Po) subspace, and the 
antisymmetric channel screens the molecular spin with a 
Kondo coupling constant 



2Vl 



|£d| + 2y|ri + 2vjr2 + s. 



(18) 



Except for the presence of molecular energy correction 
terms in the denominators, J'g and reproduce the 
Kondo couplings of Eq. ^ particularized for A = and 
/7 — > oo. The same procedure can be followed when both 
e-ph couplings are non-zero and the result of Eq. ([6]) re- 
covered. This analysis not only validates the results of 
Eq. (O and Eq. © for ujq^oji ^ Tk but it also gives some 
insight on the physics when or uji are of the order or 
smaller than Tk ■ To illustrate this we consider A = and 
uji > Tk- The dependence of Ti on 6E [see Eqs. ([T^] 
now needs to be considered and the resulting equation 
for 6E can be solved numerically. The main result is an 
increasing \5E\ (an increasing Kondo temperature) with 
decreasing uji , an effect that may be related to the high 
Kondo temperatures observed in some devices 



IV. PHONON SPECTRAL DENSITY 

Having characterized the nature of the Kondo screen- 
ing, we now analyze the molecular vibrational modes. We 
define the phonon propagators in terms of the phonon 
fields (f>Q — a + a'^ and (f>i ^ b + 

V^it) = ((0^(0,0^(0))) - -iOit) ([0^(i),0^(O)]) (19) 

for 77 = 0, 1. Calculating the self-energy to second order 
in the e-ph interaction we obtain 



2lu„ 



(20) 



where ^^{uj) indicates the Fourier transform of I?r)(i)- 
This quantity has been studied by Hewson et al?^ for a 
Holstein coupling, where the self-energy no(ijj) is given by 
the charge susceptibility. In the Kondo limit (for U > Q) 
it has little structure at low frequency and the vibrational 
mode just acquires a small width due to the coupling with 
the electrons at the molecular orbital. The self-energy for 
the CMM is given by 



ni(^) = 2VlY,mUAka,c\^J.))u 



k.cr 



+ {{c\kadcT.dlcAka))u>)- 



(21) 



Since the Kondo screening is due to the symmetric chan- 
nel, the antisymmetric electronic modes can be decoupled 
from the molecular orbital and Wiiuj) calculated in the 




FIG. 3: Phonon spectral density A{uj). a) u)\ = 0.25D, 
loi/Tk = 10, and ?/Tk = 0.15,0.25,0.35,0.45,0.55,0.6, and 
0.625 (from bottom to top in the low energy region), b) uji = 
0.025D, lui/Tk = 15, <;/Tk = 0.15,0.25,0.35,0.45, and 0.55 
(from bottom to top in the low energy region). 



lowest order in Va- 111(0;) can then be written in terms of 
one-particle propagators that for T <C Tk and lu ^ \ed\ 
are given hy -^J^k Im((cAfca, c^fc^))c^+iO = P and 



Im((d^,d;))^^+,o = t; j 

TT 1 S 



T2 ■ 



(22) 



We have 



ni(c^) 



/(^2)-/(^l) 

LO — Vl + V2 + 



bJ -\- Vi — 1^2 + iQ 

Tk/tt 



At zero temperature the integrals reduce to a simple form 
and the phonon propagator reads: 



2uji 



;2 - 2cJi^[ln(: 



- z2arctan(^)] 



with = 2V^pTK/rs = iVA/Vs)^2TK/7:. 

The phonon spectral density A{ijj) = —lmT>i(Lu)/TT is 
shown in Fig. [3] for different values of the microscopic 
parameters. For small e-ph coupling it presents a single 
peak at loi broadened by the coupling to the electronic 
degrees of freedom. For temperatures lower than Tk and 
moderate e-ph couplings, the system presents a soft mode 
(low energy peak) . For temperatures larger than Tk (not 
shown) the soft mode disappears as the low-frequency 
spectral density of the d-electrons is structureless and 
small. 

The emergence of the soft mode is associated to the 
onset of a dynamical Jahn- Teller distortion- - induced by 
Kondo correlations. To illustrate this we perform a semi- 
classical treatment in which the CMM displacement Xb 
is replaced by a real number. For simplicity we will con- 
sider only the case A = 0, but the conclusions are valid in 
the general case. The molecule is coupled to the left and 
right leads with hopping amplitudes Vl = iVo + gxb) and 
Vr = (^0 — g^b), respectively. For a given value of the 
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C-number Xb we can perform a unitary transformation 




to a new base where the {c2kcr} are decoupled from 
the molecule and play no role in the Kondo physics 
while the {cika} are coupled with a tunneling amplitude 

y/2Vo^l + -yixl that determines the Kondo screening. 

The Kondo temperature for this semiclassical treatment 
is therefore given by 

T^, ^ De-"""''''^^^'^\ (24) 

We can write the zero-temperaure energy gain of the sys- 
tem as a sum of the electronic [see Eq. (fT6|) ] and phononic 
contributions^'^i^l 

E{x') = -De-^'^P''-^^+-''^^ + ^wj^x'\ (25) 

where we have introduced the reduced displacement x' = 
^Xb- In Fig. 0] we plot the energy gain as a function 
of x' for the same parameters as in Fig. [3^. For small 
e-ph coupling, <;/Tk ^ 1, the energy has a parabolic 
form and the minimum is for x' = 0. The curvature 
of the parabola decreases with increasing and this is 
consistent with the shift to lower energies of the main 
peak in the spectral density (see Fig. [3^). For <;/Tr- > 0.3 
two symmetric minima emerge for non-zero values of x' 
showing that there is there is a Jahn- Teller instability 
associated with the formation of the Kondo effect. This 
Jahn- Teller distortion has associated a soft mode which 
is reflected in the phonon spectral density as a low energy 
peak. 




FIG. 4: Energy gain as a function of the reduced displace- 
ment x' — S-Xb of the CMM for different values of the electron 
phonon coupling <;/Tk = 0.15,0.25,0.35,0.45,0.55,0.6, and 
0.625 (from top to bottom). Other parameters as in Fig. [3^. 



V. SUMMARY AND CONCLUSIONS 



We have studied analytically electron-phonon corre- 
lation effects in molecular transistors using Schrieffer- 
Wolff transformations and a Varma-Yafet ansatz. While 
the first technique provides a description in terms of 
low energy Hamiltonians, the latter allows the study of 
the whole range of parameters, including situations in 
which the phonon frequencies are smaller than the Kondo 
temperature. We considered a Holstein coupling and a 
center-of-mass mode and showed that both produce a 
renormalization of the Kondo temperature. Remarkably, 
the Holstein coupling renormalizes the Kondo coupling 
and weakens the gate voltage dependence of Tr-. This 
may explain the puzzling behavior observed in molecular 
Co and Cu complexesi^i^ 



The center-of-mass motion of the molecule opens an 
additional channel for charge and spin fluctuations. At 
low energies the latter are described by an asymmetric 
two-channel Kondo Hamiltonian. Due to channel asym- 
metry, however, the screening of the localized spin will 
be in general dominated by a single channel. While the 
temperature dependence of the conductance may present 
some anomalies due to the additional channel, at low 
temperatures the system is a Fermi liquid and the low 
bias conductance through the molecule is G '--^ 2e^/h. 
Our calculations show that, in the Kondo regime a sin- 
gle channel is active at zero temperature and no inter- 
ference is possible between direct and phonon assisted 
channels^S. 



Finally, we showed that the phonon spectral density 
of the center-of-mass mode presents a low energy struc- 
ture induced by the coupling to the Kondo quasiparti- 
cles. Recent technical advances in single molecule phonon 
spectroscop y^^'^^ may allow for a direct study of this ef- 
fect and provide useful information on the e-ph interac- 
tion in these systems. 
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